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Abstract 

We study the effective Hamiltonian for strong-coupling lattice QCD in the case of non-zero 
baryon density. In leading order the effective Hamiltonian is a generalized antiferromagnet. For 
naive fermions, the symmetry is U{ANf ) and the spins belong to a representation that depends 
on the local baryon number. Next-nearest neighbor {nnn) terms in the Hamiltonian break the 
symmetry to U{Nf) x U{Nf). We transform the quantum problem to a Euclidean sigma model 
which we analyze in a 1/W expansion. In the vacuum sector we recover spontaneous breaking of 
chiral symmetry for the nearest-neighbor and nnn theories. For non-zero baryon density we study 
the nearest-neighbor theory only, and show that the pattern of spontaneous symmetry breaking 
depends on the baryon density. 
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I. INTRODUCTION 


The study of quantum chromodynamics at high density is almost as old as the theory 
itself [0]. In recent years the held has attracted wide interest in the wake of a revival of 
the idea of color superconductivity (CSC) |^. The stimulus for this revival was the 
observation |^, ^ that the instanton-induced quark-quark interaction can be much stronger 
than that induced by simple one-gluon exchange, and can thus give a transition temperature 
on the order of 100 MeV. Subsequent work showed that the perturbative color-magnetic 
interaction also gives rise to a strong pairing interaction. 

These and other dynamical considerations [0 underlie a picture of the ground state of 
high-density QCD in which the SU{3) gauge symmetry is spontaneously broken by a BCS- 
like condensate. The details of the breaking, which include both the Higgs (or Meissner) 
effect and the rearrangement of global symmetries and Goldstone bosons, depend on quark 
masses, chemical potentials, and temperature. Prominent in the list of possibilities are 
those of color-flavor locking in three-flavor QCD and crystalline superconductivity—with 
broken translation invariance—when there are two flavors with different densities H. For a 


review see 
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As noted, CSC at high density is so far a prediction of weak-coupling analysis. One ex¬ 
pects the coupling to become weak only at high densities, and in fact it turns out that reliable 
calculations demand extremely high densities [^. The use of weak-coupling methods to 


make predictions for moderate densities is thus not an application of QCD, but of a model 
based on it. It is imperative to conhrm these predictions by non-perturbative methods. 
Standard lattice Monte Carlo methods, unfortunately, fall afoul of well-known technical 
problems when the chemical potential is made nonzero, although we do note remarkable 
progress made recently in the small-/i regime |^, p!3[| . 

In this paper we initiate a study of high-density quark matter based on lattice QCD in 
the strong-coupling limit.^ We work in the Hamiltonian formalism, which is more amenable 
than the Euclidean formalism to strong-coupling perturbation theory and to qualitative 
study of the ensuing effective theory |^, 0.0.0- The fermion kinetic Hamiltonian is a 
perturbation that mixes the zero-flux states that are the ground-state sector of the electric 
term in the gauge Hamiltonian. In second order, it moves color-singlet fermion pairs around 
the lattice; the effective Hamiltonian for these pairs is a generalized antiferromagnet, with 
spin operators constructed of fermion bilinears. 

We depart from studies of the vacuum by allowing a background baryon density, which 
is perforce static in second order in perturbation theory. Our aim at this stage is to discover 
the ground state of the theory with this background. In third order (when W = 3) the 
baryons become dynamical; we display the effective Hamiltonian but make no attempt to 
treat it. 

The symmetry group of the effective antiferromagnet is the same as the global symmetry 
group of the original gauge theory. This depends on the formulation chosen for the lattice 
fermions. Following |^, we begin with naive, nearest-neighbor fermions, which suffer from 


species doubling [|^ and possess a global U{ANf) symmetry group that contains the ordinary 
chiral symmetries [as well as the axial H(l)] as subgroups. We subsequently break the too- 
large symmetry group with next-nearest-neighbor (nnn) couplings along the axes in the 
fermion hopping Hamiltonian. A glance at the menu of fermion formulations reveals the 


^ An early discussion of our program, with early results, was given in 0 - 
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reasons for our choice. Wilson fermions 0 have no chiral symmetry and make comparison 
of results to continuum CSC difficult if not impossible. Staggered fermions likewise 
possess only a reduced axial symmetry while suffering a reduced doubling problem. The 


overlap action is non-local in time and hence possesses no Hamiltonian; attempts 
to construct an overlap Hamiltonian directly have not borne fruit. Finally, domain-wall 
fermions ^ have been shown to lose chiral symmetry and regain doubling when 
the coupling is strong. 

As we discuss below, while the nnn theory still exhibits doubling in the free fermion 
spectrum, we are not interested in the perturbative fermion propagator but in the spectrum 
of the confining theory. We take it as a positive sign that the unbroken symmetry is now 
U{Nf) X U{Nf). This symmetry is what we want for the continuum theory, except for the 
axial U{1). The latter can still be broken by hand.^ 

Our emphasis on the global symmetries is a consequence of the fact that the gauge field is 
not present in the ground-state sector and does not reappear in strong-coupling perturbation 
theory. In other words, confinement is a kinematic feature of the theory, leaving no possibility 
of seeing the Higgs-Meissner effect directly. This is but an instance of confinement-Higgs 
duality, typical of gauge theories with matter fields in the fundamental representation . 
Our aim is thus to identify the pattern of spontaneous breaking of global symmetries. For 
various values of W and A^/, this can be compared to weak-coupling results 
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This paper is largely an exposition of formalism, along with partial results. We study the 
nearest-neighbor antiferromagnetic Hamiltonian, both with and without a uniform baryonic 
background density. We transform the quantum Hamiltonian into a path integral for a non¬ 
linear a model, where the manifold of the a field depends on the baryon background. We 
then investigate the limit of large W and show that the global U{4:Nf) symmetry is indeed 
spontaneously broken. 

Adding in the nnn couplings is a problem of vacuum alignment We do this in the 
vacuum sector and recover the result that the U{Nf) x U{Nf) chiral symmetry is broken 
to the vector U{Nf). The analysis for the finite-density theory is more involved and we defer 
it to a future publication. 

Other groups have recently studied the strong-coupling effective Hamiltonian for naive 
and Wilson fermions at non-zero chemical potential |^, |^, We differ from their ap¬ 
proaches in eschewing mean field theory in favor of the exact transformation to the a model, 
which is amenable to semiclassical treatment. As noted above, we base our program on nnn 
fermions; we also work at fixed baryon density. 

Let us walk through the paper. We review in Sec. |1| the derivation of the effective 
Hamiltonian of lattice gauge theory in strong-coupling perturbation theory [O, 17|. The 


second-order Hamiltonian [0{l/g‘^y\ is an antiferromagnet with U{4:Nf) spins; the global 
symmetry group is U{ANf) for the nearest-neighbor theory, broken to U{Nf) x U{Nf) by 
nnn terms. The baryon number at each site determines the representation of U{4Nf) carried 
by the spin at that site. In second order, baryon number is static; it becomes mobile in the 
next order, where (for W = 3) the new term in the effective Hamiltonian is a baryon hopping 
term. 

The baryon operators responsible for the hopping are composite operators that do not 
obey canonical anti-commutation relations. If this were not the case, then the effective 


^ The breaking of the naive fermions’ symmetry by longer-range terms is a feature 
0 and also occurs if naive fermions are placed on a bcc lattice [^. 


of SLAG fermions 
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Hamiltonian in third order would strongly resemble that of the t-J model 


Ht-j = cl^Cjs + J ^ 

(U> {ij) 


riiUj 


+ J'. 


( 1 . 1 ) 


Here Cjg is an annihilation operator for an electron at site j with spin s, and the number 
operators rii = cfcj and spin operators Si = ^cjcrci are constructed from it. The added 
term J'' is a more complicated hopping and interaction term. The t~J model describes a 
doped antiferromagnet; it arises as the strong-binding limit of Hubbard model, a popular 
model for itinerant magnetism and possibly for high-Tc superconductivity. The model is 
not particularly tractable and, absent new theoretical developments, does not offer much 
hope for progress in our hnite-density problem. It is nonetheless worth pondering the fact 
that a model connected, however tentatively, with superconductivity appears in a study of 
high-density nuclear matter. 

In the remainder of this paper, we work only to 0(1/5'^), where the baryons are hxed in 
position. Motivated by the similarity of our Hamiltonian to the Heisenberg antiferromagnet, 
we apply condensed matter methods developed for that problem. Indeed, condensed matter 
physicists have generalized the SU{2), spin-1/2 Heisenberg model to SU{N) in many rep¬ 
resentations ^ 1^, 1^, which corresponds to adding flavor and color 

degrees of freedom to the electrons.^ These are exactly the generalizations needed for our 
effective Hamiltonian. With W colors and N (single-component) flavors, a site of the lat¬ 
tice can be constrained to contain a color-singlet combination of mNc particles. The flavor 
indices of the spin then make up a representation of SU{N) whose Young diagram has W 
columns and m rows (see Fig. |I[). We set 


N = iNf (1.2) 

and the correspondence is complete (until we include nnn terms in the Hamiltonian). 



FIG. 1: The representation of t/(4iVj) carried by the spin in the effective antiferromagnet. m is 
related to the baryon number at the site according to m = B + 2Nf, with \B\ < 2Nf. 


In Sec. m we derive a a model representation for the partition function of the antifer¬ 
romagnet. Following Read and Sachdev [^], we employ spin coherent states to dehne 


the a held. N and m determine the target space of the a model to be the symmetric space 


^ We refer the reader to the paper by Read and Sachdev for a survey, including a phase diagram in the 
{N,Nc) plane. 
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U{N)/\U{m) X U{N — m)]; the number of colors becomes an overall coefficient of the 
action.^ As for the quantum Hamiltonian, the nearest-neighbor theory is symmetric un¬ 
der U{N) while the nnn terms break the symmetry to U{Nf) x U{Nf) (while leaving the 
manifold unchanged). 

The Nc multiplying the action invites a large-A"c analysis, and in Sec. 0 we study the 
vacuum sector, meaning zero baryon number, thereby. We return to an exercise proposed 
and solved by Smit |]^, in generalizing the vacuum sector to allow baryon number ±B on 


alternating sites; this means specifying conjugate representations of U{N) on alternating 
sites, with respectively m and N — m rows. As shown by Read and Sachdev [^, in this 
situation one can carry out an alternating U{N) rotation to convert the antiferromagnet 
into a ferromagnet with identical spins on alternating sites, and the classical {Nc = cxd ) 
analysis gives a homogeneous ground state. The result is, as one might expect, that U{N) is 
broken to U (m) x U{N — m) in the classical vacuum; the ground state energy is independent 
of m. The 1/Nc corrections to the energy do depend on m, however, and they select the 
self-conjugate m = N/2 conhguration {i.e., B = 0 everywhere) as the lowest-energy vacuum. 
Thus the true ground state breaks U{N) —> U{N/2) x U{N/2).^ When we add nnn terms 
to the action as a perturbation, we hnd that the ground state breaks U{Nf) x U{Nf) to the 
vector U{Nf), as expected. 

We turn to non-zero baryon density in Sec. 0. We study homogeneous states, in which 
all sites carry the same representation of U{N), with m > N/2. The classical vacuum of the 
a model is more elusive than for the vacuum sector, since now there are identical manifolds 
on adjacent sites but the coupling is antiferromagnetic. We begin by studying the two-site 
problem, and we learn that when one of the classical spins is hxed then when the energy 
is minimized the other spin is still free to wander a submanifold of the original symmetric 
space. If we replicate this to the inhnite lattice then we have a situation where the even 
spins, say, are hxed in direction while each odd spin wanders the submanifold, independent 
of the other odd spins. This means a ground state whose degeneracy is exponential in the 
volume, similar to some frustrated models or the antiferromagnetic Potts model |^. The 
cure to this disease comes from the 0{1/Nc) huctuations, which couple the odd spins to 
each other and make them align. In the end we hnd that the U{N) symmetry is broken by 
the vacuum to U{2m — N) x U{N — m) x U{N — m). Perturbing this ground state with the 
nnn terms is technically difficult, and we do not attempt it here despite its obvious physical 
interest. 

We close with a brief summary and discussion. The 0{1/Nc) calculation in the B ^ 0 
case is relegated to an appendix, as are other (but not all) technical details. 


^ The inverse gauge coupling 1/g^ multiplies the quantum Hamiltonian, and hence serves only to set the 
energy scale. 

® This result was obtained by Smit using a Holstein-Primakoff transformation on the quantum Hamiltonian. 
We note in passing that the 1/W calculation includes the effect of the time-derivative terms in the 
action that were dropped in the leading order, and these terms do remember the difference between the 
ferromagnet and the antiferromagnet. 
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II. THE EFFECTIVE HAMILTONIAN 


For an SU{Nc) gauge theory with Nf flavors of fermions, we write the lattice Hamiltonian 

H = He + Hu + Hf. (2.1) 

Here He is the electric term, a sum over links (n/i) of the SU{Nc) Casimir operator on each 
link, 

^ nfi 

Next is the magnetic term, a sum over plaquettes, 

1 

2c/ p 


nfi‘ 




( 2 , 2 ) 


(2.3) 


Finally we have the fermion Hamiltonian, 


'i-i 




n+/c/i,/i 




n+jfjb 


np 


i>o 


\k=0 


afS 


(2.4) 


The fermion held carries color a and flavor / at site n. The function D{j) is a kernel 
that defines the lattice fermion derivative. It can yield a naive, nearest-neighbor action if 
D{i) = a long-range SLAG derivative if D{i) = —(—or anything in between, 
such as a nnn action obtained by truncating the SLAG kernel to its j = 1,2 terms. 

For g ^ 1 the ground state of H is determined hy He alone to be any state with zero 
electric held, whatever its fermion content. 


|o)Ix)f 


_n/i 


0 ) 


\x)f- 


(2.5) 


These states have energy Cq = 0 and are degenerate with respect to all the fermionic degrees 
of freedom. We consider perturbation theory in H = Hu + He- Both Hu and He are sums 
of operators that are strictly bounded, independent of g except for the explicit coefficient 
in Hu- We can dismiss first-order perturbations by noting that Hu and He are multilinear 
in link operators U and U\ which are raising/lowering operators for the electric field; thus 
there are no non-zero matrix elements within the zero-field sector. 

We proceed to higher orders, and seek an effective Hamiltonian that acts in the zero-field 
sector |4^. Define Pq to be the projector onto the subspace of all the E = 0 states. Then 
perturbation theory in V gives an effective Hamiltonian, 

Heff = PoVQDVPo + PoVQDVQDVPo + ■ ■ ■. ( 2 . 6 ) 

Here Q = 1 — Pq projects onto the subspace orthogonal to the E = 0 states; the operator 
D = (eo — He)~^ supplies energy denominators, so that 

QD= Y. {A|. ( 2 . 7 ) 

states 
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The intermediate states |A) contain flux excitations. In second and third order the pat¬ 
terns of flux can only be strings of length j in the fundamental representation of the color 
group. Thus the energy denominators are 

^o-^x = ~9 ^Cf\j\, ( 2 . 8 ) 

where Cp = {N^ — l)/(2iVc) is the quadratic Casimir of the fundamental representation of 
SU{N,). 

The perturbations Hp and Hp are explicitly of 0{l/g^) and 0(1), respectively; each 
energy denominator gives a factor of Thus to 0{l/g^) we can forget about Hp. Our 

result to this order is 


ifeff = PoHpDHpPo + PoHpDHpDHpPo. (2.9) 

Since Hp has no non-zero matrix elements within the E = 0 sector, we have dispensed with 
Q in Eq. (|2.9|) . The first term in Eq. (|2.9|) arises for any value of Nc and is 0(1/5'^); the 
case Nc = 2 must be treated carefully, but all cases Nc > 2 are generic. The second term is 
special to = 3 and is 0(1/5'*^). 


A. Second order: the antiferromagnet 


We calculate explicitly the first term in Hes- Each term in Hp creates a string of flux of 
length j, which must be destroyed by the conjugate term. Thus 

^ 1 "), ( 2 . 10 ) 

j>0 nil PI 

where we define 

KU) = > 0 . ( 2 , 11 ) 

The matrix element of the gauge fields yields ■^SasSp'y, independent of j. 

As they appear in Eq. (p.lO|) , each is next to a "0 on a different site. This invites a 
Fierz transformation on the product of fermion fields, which we write generally as 

['iplaf.'i/ji) = . ( 2 . 12 ) 

^ A 


Here are combined site, flavor, and color indices, and we have assumed that k and 

I are always different while j and k might be equal [as in Eq. ( p.lOl) ]. The matrices T"^ 
the 16 Dirac matrices, normalized to (T"^)^ = 1 , and we have defined 


are 


4 = = ±1. 


(2.13) 


This sign factor is ±1 according to whether T"^ commutes or anticommutes with it will 
be a constant companion in our calculations. As they appear in H^g , the indices i,l are 
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the same site and color but different flavors, and likewise j, k. Leaving the flavor indices 
explicit, we obtain 


fT(2) _ 

-^eff — 


AN, 


Y, k{j)s^a 


Tl/i 


n+jfL 


- (2-14) 

j^O n 


Each fermion bilinear in parentheses is a color singlet located at a given site. The second 
term contains the baryon density® = NY'iph'if^n, and the sum J2n is the total baryon 
number B'. 

We now combine the Dirac indices with the flavor indices and write 




(2.15) 


We have dehned new matrices as direct products of the 4x4 Dirac matrices and the 
U{Nf) flavor generators, 

M^ = r^®A“, (2.16) 

and we have normalized them conventionally according to 

Ti Mm^' = ^6^^'. (2.17) 


The generate a U{N) algebra, with N = ANf. 
An alternating flip 


A. 


n* 


a, 


A 


l J 

(spin diagonalization [^) removes the matrices from the odd-j terms in Hp, 
removes the sign factors from the odd-j terms in . We have hnally 


(2.18) 
and hence 


h'A/’'VL A'(j)] B'. (2.19) 

The odd-j terms are of the form M ■ M which can be written in any basis for the U{ANf) 
algebra. The even-j terms, however, contain which is defined only in the original basis 

(CT). 


B. Single-site states 

In the zero-field sector in which we work. Gauss’ Law constrains the fermion state at 
each site to be a color singlet. The drained state |dr), with '^“f'|dr) = 0 for all (a,/), is 
the unique state with B' = 0. The other color singlet states may be generated by repeated 
application of the baryon creation operator, 

b\g... = • • •, ( 2 . 20 ) 


® This baryon number is positive semidefinite, and is zero for the drained state (see below). The conventional 
baryon number Bn is zero in the half-hlled state, and thus i?k = Bn + ‘^Nf. 
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with Nc operators tlj. (Here and henceforth, the indices combine the flavor and Dirac 

indices.) 

As noted above, at each site n the operators 

Ql = (2.21) 

generate a U{N) algebra, with N = 4:Nf. The drained state is obviously a singlet under this 
algebra. The creation operator 6^^... is in the symmetric representation of U{N) with one 
row and columns (see Fig. ^). Repeated application of h^fg... to the drained state gives 


FIG. 2; Young diagram of the representation of U{ANj) carried by the baryon operator 
the state 

lx) = 6^6^ ■ ■ ■ |dr). (2.22) 

If there are m operators 6^, then the state |x) lies in the representation with Nc columns 
and m rows (see Fig. |^). Its baryon number is B'^ = m. 

/o\ 

The second-order effective Hamiltonian preserves the baryon number on each 
site. Thus any distribution of B'^ dehnes a sector within which is to be diagonalized. In 
other words, baryons constitute a hxed background in which to study “mesonic” dynamics. 
The baryon number at each site hxes the representation of U{N) at that site, which is the 
space of states in which the charges act. 


C. Global symmetries and doubling 


The j = 
generators 


terms in Eq. (|2.19|) are of the form QnQn+p.^ 


and they commute with the 


e” = Y Ql 

n 


(2.23) 


of a global U{N) symmetry group. This symmetry is in fact familiar from the lattice 
Hamiltonian of naive, nearest-neighbor fermions: Spin diagonalization of Nf naive Dirac 
fermions transforms the Hamiltonian into that of ANf staggered fermions. In the weak 
coupling limit, there are in fact 8Nf fermion flavors—the doubling problem. This doubling 
is partially reflected in the accidental U{4Np) symmetry, which is intact in the g ^ oo limit 
and is respected by the effective Hamiltonian. Retaining terms in the fermion Hamiltonian 


/o\ 

(and thus in that involve odd separations j does not break this symmetry. 

The Nielsen-Ninomiya theorem [|^ guarantees that any fermion Hamiltonian of hnite 
range will possess the full doubling problem. This is a statement, however, about weak 
coupling only, since the dispersion relation of free fermions is irrelevant if the coupling is 
strong and the fermions are conhned. It is interesting that the accidental U{4:Nf) symmetry 
nonetheless survives into strong coupling as a vestige of doubling. 

The terms in Eq. (pi.191) with even j, on the other hand, break the U{N) symmetry, as do 
even-j terms in the original fermion Hamiltonian. It is easy to see via spin diagonalization. 
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which leaves the even-j terms unchanged, that the only generators left unbroken are the 
corresponding to 

M’' = 1 (g) A“ and 75 0 A“, (2.24) 

which form the x U{Nf)ji chiral algebra. This of course makes no difference to the 

Nielsen-Ninomiya theorem, which will enforce 8 -fold doubling in the perturbative propagator 
even without the U{4Nf) symmetry. If we are interested in the realization of the global 
symmetries of the continuum theory, though, we can study this lattice theory which has the 
same symmetry. The simplest theory one may study is thus one containing nearest-neighbor 
and nnn terms. We shall proceed to discard terms with longer range; we shall begin with 
the nearest-neighbor theory, with its accidental doubling symmetry, and later break this 
symmetry to U{Nf)L x U{Nf)fi with the nnn terms. 

Two essential differences will always remain between this lattice theory and the continuum 
theory. One is the presence of the axial U{1) symmetry on the lattice. This symmetry is 
exact, broken by no anomaly, and may make the drawing of conclusions for the continuum 
theory less than straightforward unless it is broken by hand. The other difference is the 
fact that the effective Hamiltonian for baryons (see below) is also a short-ranged hopping 
Hamiltonian. If the baryons were almost free, we would say that they are surely doubled 
like the original quarks. The fact that the simplicity of the hopping terms is only apparent, 
and that the baryons are still coupled strongly to mesonic excitations, offers the possibility 
that doubling may not return. 


D. Third order: the baryon kinetic term 

The third-order term in Hes, which only exists in the case of W = 3, is calculated via 

i/® = PqVfDVfDVfPq. (2.25) 

For a single link, we have 




and the same can be proven for a chain of links. 


1 


(')i(nc'L(nc'L(nc')7o)=s 


7^ 
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(2.26) 


(2.27) 


Thus [/, gf,... are here (temporarily) flavor indices]. 


j>0 n/i 

The kernel is 






n+jpbj yVn yVn 

imr 


+ h.c. (2.28) 


K{j) = 


. 2 • 


(2.29) 


6 (^lg^CF\j\y 

Again, spin diagonalization simplifies the odd-j terms, but not the even-j terms. The result 

nS = + hSL ( 2 . 30 ) 


IS 
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with 


(2.31) 


j>0 nfi 

j odd 

where 'r]^{n) is the usual staggered-fermion sign factor, and 

hfevL = -^ E K{j) Y, «m]/// (2-32) 

j>0 n/^ 

j even 

where Cm(^) ~ The baryon operators are 

b^ = (2.33) 

where we have written I to represent the compound index {fgh\, taking values in the 
symmetric three-index representation of U{N) (/, g,... once more combine flavor and Dirac 
indices). The odd-j part of , like that of , is symmetric under the U{N) doubling 
symmetry. The even-j part breaks U{N) to U{Nf)i x 

/o\ 

TTgg is a baryon hopping term. As mentioned in the Introduction, however, its simplicity 
is deceptive. The baryon operators 6^ are composite and hence do not obey canonical 
anticommutation relations, i.e., 


{'’t.C.VP.34) 

( 3 ) 

The separation of into a canonical kinetic energy and an interaction term is a challenge 
for the future. 


III. (T MODEL REPRESENTATION 


Because of the complexity of the third-order effective Hamiltonian, we restrict ourselves 
henceforth to the second-order theory, in which baryons are a fixed background. The theory 
dehned by HIq is a generalized spin model, with spins chosen to be in representations of 
U{N) according to the baryon distribution. We review in this section how to convert 


the spin model into a a model. The a field at each site will move in a manifold determined 
by the baryon number at that site. 


Coherent state basis 


We employ a generalization of spin coherent states [43| to derive a path integral for 


the spin model of . We recall that a given site carries generators of U{N) in a 
representation with columns and m rows, with B'^ = m. 

First we choose a basis for the Lie algebra of U{N). This consists of the generators S'*, 
with = whose matrix elements in the fundamental rep are 


( 5 ])/, = S.A 


39 - 


(3,1) 
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The corresponding charges are 


Q‘, = ’Ei’is'i’c - Inj] 

a ^ 

= ( 3 - 2 ) 


where we have subtracted a constant for convenience. The Cartan subalgebra consists of 
the operators 

= Q\. (3.3) 

We build the coherent states from the state of highest weight. The highest-weight state 
I To) in the representation is an eigenstate of the Cartan generators, 




(W/2)|To) fori = l,...,m 
— (W/2)|To) for i = m + 1 ,..., W 


(3.4) 


In this state, the generators take the simple form 

(ToigJiTo) = (3.5) 

with 

The state |To) is invariant (up to a phase) under the subgroup of U{N) that commutes with 
A; this is U{m) x U{N — m). The most general rotation of |To) is carried out with the 
generators Q* that are not in the corresponding subalgebra. We choose parameters a^, with 
A G [1, m] and p G [m + 1, A^], and write 


I a) = exp 



Qx - I I'l'o)- 


(3.7) 


The only generators that appear in Eq. (|3.7| ) are those that lower an Hi that starts from 


Nf./2 in Eq. (|3.4|) while raising another that starts from Any other generator 

would annihilate |To) and thus give no effect in the exponential. 

The coherent states are normalized. 



(a a) = 1, 

(3.8) 

and over-complete. 

J da a)(a = 1. 

(3.9) 

In Eq. (|3.9|) the integral is over the coset space 17(iV)/[t/(m) x U{N — m)] 

(see below). 

Matrix elements of the generators 

are given by 



(fllQjlci) = —NcCij, 

(3.10) 

where the matrix aij is given by a 

unitary rotation from A, 



a = U {a)AU (a)^ 

(3.11) 
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The matrix U{a) is built out of the m x (A^ — m) matrix a^, 



/ 0 a\ 

exp 

-1 

0 

1 

_ 1 


(3.12) 


a is both Hermitian and unitary. 

The manifold of matrices a is the coset space U{N)/\U{m) x f/(iV—m)], a sub-manifold of 
U{N). This is because for any matrix 17(a), one can generate an orbit U{a)V by multiplying 
with a matrix 

(3.13) 

where X G U{m) and Y G U{N — m). All matrices in the orbit will give the same matrix 
a when inserted into Eq. (|3.11| ), and thns in integrating over the conhgnration space of a 
one mnst choose only a single representative of each orbit. This set of representatives, the 
coset space U{N)/\U{m) x U{N — m)], is the qnotient space of the non-invariant snbgronp 
U (m) X U{N — m). 

The measnre over the coset space mnst be invariant nnder nnitary rotations. 


a) R{V)\a), 


(3.14) 


where 72(17) represents the rotation 17 in Hilbert space. Eqnation (|3.1(]|) shows that this is 
a rotation 


a —>• V aV^ 


(3.15) 


and by Eq. (|3.11| ), this means that a measnre in U mnst be invariant nnder U —> VU. This 
Exes the measnre nniqnely to be the Haar measnre in U{N), and thns one can integrate over 
the coset space by integrating with respect to U over U{N) and using Eq. ( |3.11|) . 

A representation whose Yonng diagram has N — m rows is the conjngate to the repre¬ 
sentation with m rows. Its coherent state space can be constrncted to look the same, with 
only a sign difference. To do this we start with the /owesf-weight state, which satishes 
[cf. Eq. (U)] 


77,1 To) 


-(Yc/2)|To) for2 = l,...,m 
(Ai'c/2)|To) for 2 = m-I-1,..., Y. 


(3.16) 


This introdnces a minns sign into Eq. 
the replacement of Eq. (^.10) by 


). The snbseqnent steps are identical, with only 


(a|Q)|a) 


2 ■ 


(3.17) 


Here, too, a is given in terms of A and U by Eq. ( p.ll| ). 


B. Partition function 

The partition fnnction Z = Tre“^^ can be written as a path integral by inserting the 
completeness relation (|3.9|) at every slice of imaginary time. This gives 

Z = y Da{T) exp—S', (3.18) 
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where the action is 





1 — (a(r)|a(r + dr)) 
dr 


+ H(a(T)) 


(3.19) 


The Hamiltonian H{a) is a transcription of the quantum Hamiltonian to the classical a 
matrices. Starting with the quantum operator we have 

Ql = 

= M^.gjH + ^iVeTrMT (3.20) 

Expressed in these variables, the quantum Hamiltonian is^ 

h'S = Y. JiQlQYif. 

nfj, 

9#0 


where Jj = {2/Nc)K{j). We transcribe this according to Eq. ( p.l0|) to obtain the classical 
Hamiltonian, 

. , -r 9 


H{a) = 


f) 

" / nfi 

jVO 






where 


a. 


1 = TiM^'^ar 


(3.22) 


(3.23) 


Recall that each Uri is an x iV matrix ranging over the coset space appropriate to site n. 
The time-derivative term in S' is a Berry phase [^. It can be expressed in terms of the 


matrix U that determines a via Eq. (|3.11|) . The result is^ 


S = 


dr 




(3.24) 


will vary from site to site if m does. If one takes the route of Eq. (|3.17|) for a representation 
with N — m rows, then the kinetic term for that site acquires a minus sign (see below). 

The number of colors has largely dropped out of the problem, since a is just an N x N 
unitary matrix held. The explicit factors of W in Eqs. (p.22 ) and ( p.24| ) invite a semi classical 
approximation in the large-W limit. This of course neglects the A"c"dependence of the 
couplings Jj ~ l/A)?, but the common scale of the couplings only serves to set an energy 
scale. The ground state will be independent of this scale, although correlations and the 
temperature scale will rehect it. We take the point of view that after all W = 3, and we are 
interested in properties of the effective theory for this value only. The ^ oo limit, for 
hxed couplings Jj, will be a device for investigating the properties of a generalized effective 
theory. 


7 

8 


The B' term from Eq. ( ^.19| ) indeed disappears. We have dropped an additive constant that is independent 
of B'. 

This is correct only if U is of the form given in Eq. ( 3.12| ), and in that case U cannot be integrated over 
all of U{N). 
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IV. ZERO DENSITY 


The simplest way to set the baryon density to zero is just to choose Bn = 0 on each site, 
meaning m = N/2 = 2Nf. It turns out to be just as easy to consider a slightly generalized 
case 1^, in which Bn is chosen to alternate, Bn = ±&, on even and odd sublattices. This 


means to choose a representation with m = N/2 + b rows on even sites and N — m = N/2 — h 
rows on odd sites, which gives a pair of conjugate representations of U{N). In view of 
Eq. ( p.lTp , we can substitute a —a on the odd sites and thus have identical manifolds on 
all sites. The Hamiltonian is then 


= T.Ji~ 


■lyaio" 


Tl/i 

3^0 


n+jfj. 


(<) 


i+1 


(4.1) 


m is a new parameter in the theory, and we can ask what value of m gives the lowest 
energy for the ground state. We will see that m = N/2 is indeed preferred, but only in the 
next-to-leading order in 1/iVc. 


A. Large-iVc limit 

In the large-W limit, we seek the classical saddle point® of S. We assume the saddle is 
at a conhguration (T„(r) that is independent of time, and so we drop the time derivative. 
We begin with the nearest-neighbor Hamiltonian,^® 

= (4.2) 

Again, the matrices an are Hermitian and unitary, and the expansion coefficients a/^ are 
real and satisfy = N/2. The minimum of H is clearly at a constant field, an = o'q, 

which can be diagonalized to cJn = A by a U{N) rotation. This is a Neel state in the original 
variables. The U{N) symmetry is broken to U{m) x U{N — m) and there are 2m{N — m) 
Goldstone bosons. 

The classical energy density (per link) is Cq = —JN/2, independent of m. Thus at leading 
order in 1/W, the optimal value of m is undetermined, and any alternating background of 
baryon number is equally good. 

Since it turns out that the l/W corrections select m = N/2, let us consider the effect of 
the nnn term in H for this case only. The perturbation is 

(4-3) 

n/j 

It breaks the U{N) symmetry to SU^Nf)^ x SU{Nf)ii x U{1)a x U{1)b- Assuming that 
J' -C J, we again seek the minimum energy configuration in the form of a constant field; we 
minimize 

e' = ^ (4.4) 

firi 


9 

10 


Note that the kinetic term is pure imaginary. 

This classical Hamiltonian is ferromagnetic; the antiferromagnetic nature of the quantum Hamiltonian 


(3.21) is preserved by the alternating signs in the time-derivative term. 
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among the [/(A^)-equivalent cr = cxo states that minimize the nearest-neighbor action. It 
is not hard to show (see Appendix 0 ) that e' is minimized for a = 70 0 1 . This is a 
condensate that is symmetric under the vector generators = 1 ® A“ but not under the 
axial generators 75 ( 8 ) A“, and thus it breaks the chiral symmetry to the vector subgroup, 
SU{Nf)v X U{1)b. 


B. 1/A^c corrections 


Returning to the nearest-neighbor theory, we consider fluctuations around the an = A 
minimum of S. First we rescale r — 2 r/Ac in order to put the kinetic and potential terms 
on an equal footing, giving 


S 


2 



n ^ rifi 


( 4 . 5 ) 


with [3 = Ncj3/2. Recalling Eq. (|3.12|) , we can write 




( 4 . 6 ) 


where An is anti-Hermitian and anticommutes with A. It is more convenient to work with 
the Hermitian matrix 


Ln= 2A„A, 


( 4 . 7 ) 


in terms of which we expand 

an = A -|- Ln — —L^A. (4-8) 

If we further expand in the basis of generators of U{N), 


Ln = Y.^'^M\ 

V 


( 4 . 9 ) 


we hnd that the corresponding to generators of U{m) x U{N — m) vanish; this is the 
subgroup under which the vacuum is symmetric. The held Ln thus contains 2m{N — m) 
real degrees of freedom, corresponding to the Goldstone bosons. 

We expand Un and an in powers of Ln] using Eq. (|4.9|), we obtain to second order 


S = So + 


Nr 




1 W 7 

O <'n^T<'n ' o 2^ {^n+n L 

8 n 


8 


(4.10) 


The coefficient matrix is 

C-w'= Tr (A[M^M'^']), (4.11) 

and the classical energy is S'o = — 7 NNsd(3. C is antisymmetric and purely imaginary; 

we show in Appendix ^ that C has eigenvalues ±1, each with degeneracy m{N — m). We 
change basis so as to diagonalize (7, and write the index rj as the compound (a, ±) with 
a = 1,..., m{N — m) and the ± corresponding to the eigenvalue of C. Since the original 
are real, we have 

(/“+)* = /“-. (4.12) 
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Thus we eliminate /" and write 


Nr 


S = S„ + ^/<iri: 


-inimC’S.C + ^IZIC 


ja+ 

n+fi 


(4.13) 


The alternating sign in Eq. (|4.13|) is what makes the theory antiferromagnetic. It forces 
us to differentiate between even and odd sites, and we transform to momentum space as 
follows (dropping the + superscript): 


= 



E 


^ even, 


(4.14) 


The even sites comprise an fee lattice with lattice constant 2, and the momenta k take values 
in its Brillouin zone. We obtain 


N, 


s = So+ ^J: drill i;)%M{k) 


h 

I, 


(4.16) 


Here 


M{k) = 


^ Jd — dr —Jd'y{k) \ 

^ —Jd'){k) Jd + dr J 


and 7 (fc) = ^E/^cosfc/,. 

The gaussian path integral over the I held now gives the free energy. 


Nr 


F = Fq + m{N — m )—^ ^ 


(5 2 2 


(4.16) 


(4.17) 


where uj{k) = dJNc\Jl — j'^{k) and Fq = —\JdNNs{Nc/2Y. For the ground state energy, 
we take /3 —> cx) and obtain (restoring all constants) 


Fa = -JN,Nd 


Nr 


1 + 


1 m{N — m) 


Nr 


N 


Ibz I 2tt 


1 - Jl--f^{k) 


(4.18) 


This is exactly the result of Smit |^. The 0(l/iVc) corrections lift the degeneracy of the 
ground states with different values of m. The integral in Eq. ( [4.18|) is positive and its 
coefficient contains the number of Goldstone bosons. Thus the state of lowest energy is that 
with m = N/2, and the symmetry breaking scheme is U{N) —>■ U{N/2) x U{N/2). Further 
breaking by the nnn terms was discussed above. 


V. NONZERO BARYON DENSITY 

The zero-density theories considered in the preceding section were dehned by selecting 
representations with m and N — m rows on alternating sites. For any m, this led to a cr 
model with identical degrees of freedom on all sites—after redehnition of the spins on the 
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FIG. 3: Two U{N) spins in different representations 

odd sublattice—and ferromagnetic couplings.We eventually settled on m = N/2 as the 
background that gives the ground state of lowest energy. 

Introducing non-zero baryon density means changing m on some sites of the lattice. Since 
in general there will be representations on different sites that are not mutually conjugate, 
different sites will carry a variables that do not live in the same submanifold of U{N). We 
here limit ourselves to the simpler case of uniform m, where adjacent sites carry identical 
spins—but the coupling is antiferromagnetic. 

In order to learn how to work with such a theory, we begin by studying the two-site 
problem. The results of this study will lead directly to an ansatz for the ground state of a 
lattice with a hxed density of baryons. 


A. The two-site problem 

1. Classical solution 


Consider, therefore, two sites with quantum spins Qi and Q 2 that carry representations 
of U{N) with mi and m 2 rows, and W columns (see Fig. The quantum Hamiltonian is 

H = JQIQI (5.1) 

an antiferromagnetic coupling. The corresponding classical a model has the interaction 
Hamiltonian 

,7 / N \ 

^ 2\y) (5.2) 

where 

(Tj = UiKiU}. (5.3) 

The two A matrices reflect the different values of m* according to 


w = 



0 

1 V—mi 


(5.4) 


11 


The (—1)" factor in the kinetic energy retained information about the antiferromagnetic nature of the 


quantum problem; it did not affect the classical analysis. 
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The Nc ^ oo limit is the classical limit, in which we seek values of ai ^2 that minimize 
H{a). A global U{N) rotation, viz., 


a, ^ VW*. 


( 5 . 5 ) 


can be used to diagonalize ai so that cti = Ai. Now we have to minimize Tr A 1 CT 2 . The case 
of conjugate representations, m 2 = N — mi, is easy: a 2 is a unitary rotation of A 2 , which 
(in this case) can be rotated into — Ai. This is the unique antiferromagnetic ground state. 
cTi and (72 can be copied to the odd and even sublattices of an inhnite lattice to give the 
classical Neel state considered in the preceding section. 

The case mi = m 2 = m is more complex. We consider m > N/2 for definiteness; the 
other case is similarly handled. We write explicitly^^ [from Eqs. (|3.11|)-(|3.12|)] 


0-2 




( 5 . 6 ) 


a^a is a square matrix of dimension N — m and aa^ is a square matrix of dimension m. Since 
a 2 is a rotation of A, 


2 m — A^ = Tr (72 



and hence the energy is 


E 


J 

2 

J 



Tr A(7 
Tr cos 



+ 2m — N 


( 5 . 7 ) 


( 5 . 8 ) 

( 5 . 9 ) 


E is minimized when all the eigenvalues of a^a are equal to 7r^/4. This means that the N — m 
column vectors a* form an orthogonal set in m dimensions, with 


a 


a, = 



(5.10) 


Since m > A^ — m by assumption, such a set of vectors can always be found. 

Since a^a = (7r^/4)liv-m, we have sin (2\/ ala^ = 0 and so the off-diagonal blocks of 
Eq. (|5.6|) vanish. The lower-right block of (72 is the unit matrix Iw-m- We know that 
(<^ 2 )^ = 1 since A^ = 1 and thus the upper-left block must have eigenvalues ±1. Equating 
traces of (72 and A, we find that the upper-left m x m block must take the form 

where 

A(”*) = f ^2m-iV 0 \ 

\ U —iN-mJ 


This generalizes a parametrization found in . 
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and G U{m). represents the coset U{m)/\U{2m — N) x U{N — m)]. 

We conclude that the classical ground state of this B ^ 0 two-site problem is degenerate, 
even beyond breaking the overall U{N) symmetry. The solutions can be written as 


0-1 



0-2 


\ 0 iN-m ) 


(5.13) 


(to which a global U{N) rotation can be applied). A particular instance of is 
given by Eq. (|5.12|) . The symmetry of the vacuum is the set of rotations that leaves both 
(Ji and (72 invariant, namely, U{2m — N) x U{N — m) x U{N — m). 


2. Quantum fluctuations 


The classical solution of the two-site problem will guide us in approaching the problem 
of an inhnite lattice below. We expect that spontaneous symmetry breaking will give a 
vacuum of the same character, with continuous degeneracy. There are, however, two kinds 
of degeneracy in the two-site problem: that which results from breaking the global U{N) to 
U{m) X U{N — m), and that which comes of breaking the U{m) subgroup to U{2m — N) x 
U(N — m). The latter degeneracy is connected with freedom in choosing the orientation of <72 
relative to ai. It is instructive to see how quantum mechanical fluctuations lift the classical 
degeneracies. 


The quantum two-site problem is easy to solve. We rewrite the Hamiltonian (5.1) as 


h = 2 

2 L 


(Qi + Q2Y ~ Qi ~ Q 


(5.14) 


and Q 2 are constants, the quadratic Casimir operator in the m-row, W-column repre¬ 
sentation of U{N). The first term in Eq. (^.14|) is minimized by coupling Qi and Q 2 to 
the representation that minimizes the Casimir, which is the representation with 2m — N 
rows and W columns (see Fig. |^). The ground state has discrete degeneracy equal to the 
dimension of this representation. 

The exact quantum solution naturally shows no sign of spontaneous symmetry breaking 
and hence it is not of much relevance to the infinite volume problem. More interesting is the 
problem where the state of Qi is fixed and Q 2 is allowed to vary. This breaks by hand the 
global U{N) while allowing quantum fluctuations to lift any remaining degeneracy in the 
relative orientation of the two spins, so it can be regarded as quantization in the presence 
of spontaneous symmetry breaking. In effect, this is mean held theory. 

We replace the Hamiltonian (|5.1|) by 


j^MF 






(5.15) 


To minimize the energy we maximize the mean held by choosing the state of Qi to be the 
highest-weight state. This state diagonalizes the generators Hi of the Cartan subalgebra 
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FIG. 4: Coupling two spins in a 7-row representation of 17(12) to the representation with minimal 
Casimir operator 


while other generators of U{N) have expectation value zero. Thus 

N 

= JY.{Hu)H2,. (5.16) 

i=l 

The operators all commute, and their eigenvalues make up the weight diagram of the 
representation.^^ jg p^-oduct of the weight vectors of the two spins. The energy 

is minimized by choosing for Q 2 a state that lies opposite the highest-weight state in the 
weight diagram. As shown in the example of Fig. ^ this still leaves a degeneracy, albeit a 
discrete one. We stress that this degeneracy comes from freedom in the relative orientation 
of Qi and (^ 2 ; it remains after quantum fluctuations lift the continuous degeneracy of the 
classical system. 



FIG. 5: Weight diagram for A = 3, m = 1, Ac = 2 [the sextet of 517(3)]. The highest-weight state 
lies at the top of the triangle; all states at the base minimize . All values of Ac give the same 
triangular shape, with multiplicity of one along the boundary. There are Ac -|- 1 states at the base. 
For A = 4 the triangle becomes a tetrahedron, and there are (Ac -|- l)(Ac -|- 2)/2 states at its base. 


More precisely, the weight diagram shows eigenvalues of the A — 1 traceless diagonal generators of SU (A). 
These can be obtained by isolating the 17(1) member of the set Hi and taking linear combinations of the 
rest. 
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In the iVc —>■ oo limit, the discrete degeneracy becomes inhnite and presumably it is well 
described by the continuous degeneracy of the classical problem. 


B. Infinite lattice 


At A'",; = oo we seek the saddle point of the action, which we assume to be a time- 
independent conhguration. The classical Hamiltonian of the a model is 

H = (5.17) 

Seeking an antiferromagnetic ground state, we set cr^ = A on the sublattice of even sites. 
The odd sites are then governed by 

_ Jd ^ TrAcx^. (5.18) 

n odd 


This is just the two-site problem studied above, replicated over the lattice. As we saw above, 
the ground state conhguration is degenerate with respect to the conhguration at each odd 
site, 

° ' (5.19) 


^n. 



1 V-' 


N—nri 


A uniform choice for the odd sites, = A*-"*^ for instance, breaks the U(N) symmetry to 
U{2m — N) X U{N — m) x U{N — m); a non-uniform choice can break the symmetry all 
the way to U{N — m). The entropy of this classical ground state is evidently proportional 
to the volume. 

As noted for the two-site problem, the continuous degeneracy of the ground state is an 
artifact of the classical, Nc ^ oo limit. Quantum huctuations will spread each odd spin’s 
wave function over the U{m)/[U{2m — N) x U{N — m)] manifold. A mean-held ansatz for 
the even spins, as noted in the discussion of the two-spin problem, will still leave a discrete 
degeneracy for the odd spins; the symmetry breaking scheme will depend on how the odd 
spins are allocated to the available states. Furthermore, one can contemplate making a 
non-uniform ansatz for the even spins as well, reducing the phase space available for odd 
spins with unequal neighbors but adding entropy on the even sublattice. The situation is 
reminiscent of that in the antiferromagnetic Potts model [^, the phase structure of which 
is not yet understood. 

This Potts-like discrete degeneracy, however, is an artifact of the mean-held approach 
that, like the classical ansatz, assumes a hxed state for the spins on the even sites. An 
essential difference between our a model and the Potts model is that our degrees of freedom 
are continuous and will huctuate as soon as they are allowed to do so. A given odd spin 
will not be surrounded by a uniform hxed background; the neighboring even spins will be 
inhuenced by all their odd neighbors, and will induce an interaction among the odd spins 
that makes them rotate together. This should reduce the entropy of the ground state to zero. 
The systematic way to see this ehect is to carry out a l/A). expansion around the classical 
ansatz, which we do in Appendix p. The result is a ferromagnetic interaction among the 
an on the odd sites. Thus the ground state turns out to be the two-site solution, replicated 
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uniformly over the lattice: 




A, 

/ 0 \ 

Iv 0 J ’ 


n even, 
n odd, 


(5.20) 


where is a global degree of freedom. The symmetry group of the 

vacuum is U{2m — N) x U{N — m) x U{N — m). 


VI. SUMMARY AND DISCUSSION 


Let us summarize the results presented in this paper. In the vacuum sector, we have 
rederived Smit’s result for the lowest-energy conhguration of alternating Bn = ±(m — N/2) 
sites. The result is indeed Bn = 0; the U{4Nf) symmetry of the nearest-neighbor theory is 
spontaneously broken to U{2Nf) x U{2Nf). We extended this result to the nnn theory and 
found that its U{Nf) xU{Nf) chiral symmetry is broken to the vector U{Nf) flavor subgroup. 
Adding net baryon number to the system, we examined the case of uniform baryon density. 
Bn = m — 2Nf > 0. Here our study was limited to the nearest-neighbor theory, and we 
found a Neel-like ground state that breaks U{N) to U{2m — N) x U{N — m) x U{N — m). 
The number of Goldstone bosons ugb thus depends on the baryon density Bn as 


^gb — 2{2Nf — \Bn\){2Nf + 3|i?„|). 


( 6 . 1 ) 


Directions for future work begin with adding nnn interactions to the B > 0 theory and 
extracting from it a prediction for the breaking of the continuum-like chiral symmetry. (The 
axial U{1) symmetry can be broken by hand.) Another direction is to gain greater freedom 
in hxing the baryon density. A constant value of Bn > 0 means a baryon density that is close 
to the maximum allowed on the lattice; the density can be lowered by setting Bn ^ 0 only 
on a sparse sublattice, along the lines shown in [|^. An ultimate goal, as for the Hubbard 
model, is the incorporation of the third-order term in the effective Hamiltonian in order to 
have a theory with dynamical baryons. Perhaps an instructive half-measure would be to 
study the second-order theory in the presence of a disordered baryon background. 

The strong-coupling effective theory can be regarded as a QCD-like model, possessing 
gauge invariance and the correct degrees of freedom. In that case the lattice spacing is 
merely a parameter, an overall scale. More insight can be gained by considering the strong¬ 
coupling theory to represent QCD at large distances, derived by some renormalization-group 
transformation from a weak-coupling short-distance Hamiltonian. On the one hand, one 
would expect any such effective Hamiltonian to contain many terms of great complexity; on 
the other hand, a simple lattice theory such as ours might offer a qualitative approximation 
to the real theory (as long as one accepts the loss of Lorentz invariance). We can estimate 
the lattice spacing to be some scale at which the running QCD coupling is large, certainly 
greater than the radius of a proton. The limitation that the lattice puts on the density 
then becomes a physical issue. Taking the lattice spacing a to be on the order of 1 fm, 
the highest baryon density allowed by the lattice is 2Nf fm“^. For low values oi Nf this 
may not be enough to see hnite-density phase transitions, in particular a transition to color 
superconductivity. Perhaps a way out is to consider an unphysically large number of flavors. 
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APPENDIX A: MINIMIZING THE NEXT-NEAREST-NEIGHBOR TERM 


P 


The signs are dehned only when the are written in the basis 

A“. We choose a chiral basis for the gamma matrices, so that 75 = 


cr^ 


and /3 = ph The energy is a sum of squares. 


= P^ ® = 


V 


(Al) 


with the constraint = N/2. The coefficients take on the values 

{—3,—1,1,3}. The minimum of e' occurs when all are zero except those correspond¬ 
ing to Ar, = —3, namely, those for = /I (g) A“ = (g) A“ and = (3'y^ (g) A“ = p^ (g) A“; 

the energy is independent of these Thus the set of solutions can be written in the form 


o-Q 


/ 0 0 R 0 \ 

0 0 0 H 

Rt 0 0 0 

V 0 f/t 0 0 y 


^ TT , P^ -ip^ ^ rrt 
-Z— ®u^ --- ®U\ 


(A2) 


Recalling that cxg = 1, we have HRt 
with 

/f/ 0 0 0 \ 
0 R 0 0 
^ “ 0 0 10 
Vo 0 0 ly 


1, so U G U{Nf). A chiral rotation cxo — V^aoV, 


^(l + p3)®R+l(l + p3)0l, 


(A3) 


turns (To into p^ = Jo, which is invariant only under vector transformations generated by 
1 0 A“. 


APPENDIX B: THE MATRICES C AND 


To analyze the matrices C and Dn, given by Eqs. ([4.11|) and ( |C18|) , we begin with the 
U{N) generators M"^ that lie outside the subalgebra H = U{m) x U{N — m) that commutes 
with 


A 


0 


-1 


N—m 


(Bl) 
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We choose for them a basis and with p = 1,... ,m and g = m + 1,..., iV, given 
by 


(A^’’»‘)/« = i (KfS„ + SM (B2) 

^ (SvA, - S,Al) ■ (B3) 

Since the coset space U{N)/H is a symmetric space, the commntator lies in 

H; in order for Tr '^ “ ]) to be nonzero, the commntator mnst have a nonzero 

component in the Cartan snbalgebra of H. This is only possible if a 7 ^ a' and (p, q) = (p', q'). 
Thus in this basis C takes the form 

7-y • I 0 lm(iV—m) 

\ — Im(V-m) 0 

Diagonalizing C gives eigenvalues ±1. The generators corresponding to the basis that diag¬ 
onalizes C are 



= S„S„ (B5) 

{M”-)s, = (M«‘ - UPAfs = S^Ar (B 6 ) 

As for D: The only anticommutators among the that do not vanish [note the 

bounds on (p, q)] are between and ~, viz., 


MP }fg - SqqiSpf6p>g + Spp'SqffSqg. 

Noting that we hnd that Dn takes the block-diagonal form 


jjpq±,p'q'± _ 


E TV 

m{n) 

= V E (iv 

m{n) 


_ 


a. 


(m) 

m 

0 -lv_ 


N—m 


m JPP 


We summarize this by writing 

D = ( ® 

where En is an m x m matrix given by 

= E ( 1 - - ^ 

m(n) 


(B7) 

(B 8 ) 

(B9) 

(BIO) 

(Bll) 


It is easy to prove that the eigenvalues of E range from 0 to Ad. In particular, D is positive. 
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APPENDIX C: l/N^ CORRECTIONS TO THE R / 0 PROBLEM 


We build on the Nc = oo vacua described in Sec. 0 by allowing fluctuations around 
them. We let the cr„ on the even sites fluctuate around A; we let the on odd sites roll 
freely around the U{m)/[U{2m — N) x U{N — m)] manifold covered by Eq. (|5.19|) , and also 
execute small oscillations off the manifold into the U{N)/[U{m) x U{N — m)] coset space. 
Our goal is an effective action for the classical part ( b.l9| ) of the odd spins. To reach this, we 
integrate out the even spins; the off-manifold fluctuations of the odd spins must be included 
for consistency in the 1/W expansion. 

The counterpart of the action ( |4.5| ) for our problem has an antiferromagnetic spin-spin 
interaction, with no (—1)”' factors. We separate it into odd, even, and coupled terms, 

S = Y > 

= /dr ^ ^IAUl^rUr^, 

n odd 

= j dT Y. TrA[/ta,t/„, 

n even 

s-’" ^ jdT Y. (Cl) 


Here cr° = J2m(n) where m{n) are the nearest neighbors of the even site n. We expand 
the held on the even sites around cxn = A in the manner of Eq. (|4.8|) , 


an = A + Ln — 2 ^^^ even), 

while for the odd sites we write (see Appendix 0) 


(C2) 


ar 


= u, 



-A' 


V + \L'‘^A' 


Ul (n odd). 


(C3) 


with 


and 



Un 


f/M 0 \ 
0 Iv —m j 


(C4) 

(C5) 


Ln describes the fluctuations of the even spins around their classical value A. Un rotates the 
odd spins within the manifold of their classical values, while L'^ describes their fluctuations 
outside that manifold. We further define 



0 \ 

0 lAT-m j ’ 


(C6) 


the classical held on the odd sites. 
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We leave alone and expand and around the classical values of the helds, 


= —JdT Y. TiM^drL^, 


(C7) 


= So + 


J 


dr Y (Tr - xTr Ll^Aan - Tr L„L^ 


odd 


+dJ I dr (-TrAi„ + -TrL; 


(C8) 


Here 

in = 1 ° ) Ui (C9) 

is the rotated fluctuation held on the odd sites, and the Hermitian matrices dn and are 
sums over the odd neighbors of the even site n, 

ffn = E (CIO) 

m(n) 

Lr.= Y Lrr^■ (dl) 

m(n) 


Since both ct„ and A are block diagonal, the hrst trace in each integral in Eq. 
Now we organize the partition function as follows: 


IS zero. 


Z = 


dan exp 


odd 


A(s°" + S„+(^/dr E Tr£ 

^ \ ^ n odd 


(C12) 


with 


7 _ 

•^even 


Y[ din] exp 


Y 
' 2 


' dr 


X Y ( “iTr ALndrLn - jTr L^Acr^ - ^Tr 


(C13) 


Zeven IS a product of decoupled single-site integrals. Again we expand in the group algebra. 


Lr. = IIM\ 


(C14) 


where the sum is over the 2m{N — m) generators of U{N) that are not in U{m) x U{N — m). 
Ln can be expanded similarly and we obtain the following form for the integral over the 
even helds: 


Din exp 


Y 

2 



(CIS) 


where 


8 8 


The matrix C is the same as in Eq. (|4.11|), 


(C16) 




(C17) 


27 











while the new matrix D varies with the site n according to the average of its neighboring 
spins, 

dW = _Tr {{M\ (C18) 

We study the two matrices in Appendix O. Diagonalizing C as before, we arrive at 


nq 


^even = H W ~ 




dr 




(C19) 


where Adrj is the m x m matrix 


Xin — ^ {dr + JEn) 


(C20) 


The quantities (and E) for each q = 1,..., N — m are complex m-component vectors; they 
are rotations of the 2m{N — m) real components E (and E) into the basis that diagonalizes 
C. The matrix En is given by Eq. ( |B11|) ; it carries the dependence on an- Performing the 
gaussian integration we get 


Egven 


nq 


Det Ad, 


exp 


2 


9' 


drWMXll 


(C21) 


Finally we separate the integral (|C12|) over the odd spins into an integral over the classical 
held a^ and an integral over the huctuations around it. We obtain 




Z= / exp-^ I S'o + (A^ - m)^Tr logAd, 


X / Dim exp - 




+ / dr 


^ mq \ ^ / jiq 


(C22) 


Here m stands for an odd site, n for an even one. 

Equation (|C22|) gives an effective action for the classical odd spins a^- These enter the 
exponents through jCi [via Eqs. ([Bll|) and (|C2U|) ] and through [via Eqs. (|C5|) and (|(J11|) ]. 


The action in the hrst exponent is minimized when each matrix Enic^m) the largest 
number of zero eigenvalues, each of which makes Tr log Ad,, approach —cx). It is easy to 
check that En has 2m — N zero eigenvalues (the maximal number) when the al^ on o-d the 
odd sites m{ri) align, i.e., 


0-^ 


= do e U{m)/[U{2m - N) x U{N - m)]. 


(C23) 


Moreover, when Eq. ( C25 ) holds, all the /^’s align parallel to each other and is maximized; 
also the eigenvalues of Ad“^ are maximized (to +oo). Thus the action in the second exponent 
also has a minimum at this point in conhguration space. These effects add up to an effective 
ferromagnetic interaction among the 2d nearest neighbors m of each even site n. This 
effective interaction will align the classical spins on the odd sublattice to the same direction 
in their submanifold, U{m)/[U{2m — N) x U{N — m)]. 

The divergences in the effective action have their origin in the fact that the semiclassical 
corrections are calculated as gaussian integrals in the even fluctuation helds and the 
coefficient matrix Ad„ acquires zero eigenvalues. The correct range of integration over 
is of course not inhnite, but rather the volume of the U{N)/[U{m) x U{N — m)] manifold. 
This will regulate the divergences, but leave the effective action for the odd spins attractive. 




















APPENDIX D: FLUCTUATIONS ON THE ODD SITES 


In the classical analysis, the helds on the odd sites take values in the sub-manifold 
U{m)/\U{2m — N) x U{N — m)] of the manifold U{N)/\U{m) x U{N — m)]. We denote 
these values 


\ 0 iN-m ) 


(Dl) 


Here 


^{m) _ jj{m) 


with G U{m), and 




^2m—N 

0 Iv—m 


(D2) 

(D3) 


contains 2[N—m){2m—N) independent degrees of freedom. Any a G U{N)/\U{m) x 
U{N — m)] can be written as 


a = 


V 


cos ( 2 Vaa^'j 
sin (2 V at a) 

V a^a 


in ( 2 V aia^ ^ 
V a^a 
)s (23 


(D4) 


[cf. Eq. (|5.6|) ], which coincides with Eq. (|D1|) if 


= 


0 


{TT/2)lN-r 


(D5) 


Recall that a is an mx [N—m) matrix, so the zero block has dimensions {2m — N) x {N — m). 
We allow motion out of the sub-manifold by allowing a to vary further. 


_ jj{m) 


0 


(D6) 


The 2{N—m)^ degrees of freedom in d complement the 2{N—m){2m — N) degrees of freedom 
inherent in to give 2m{N—m), the dimensionality of the entire U {N) / [U {m) xU {N—m)] 
coset space. 

Writing a with the generalized a we have 


a = U 


( l2m-7V 
0 


V 


0 


with 



UK 


(D7) 


(D8) 
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We can also write this as 


a = U 


l2m-Af 0 

0 


[/t. 


(D9) 


fj[ 2 {N m)| jg ^ matrix in the manifold U{2{N — m))/\U{N — m) x U{N — m)]. Indeed for 
a = (7r/2) In- mi we have 


^{2{N-m)] ^ f Iw-m 0 \ ^ _jy 

y U i-N—m J 

Since U{2{N — m))/[U{N — m)x U{N — m)] is a self-conjngate manifold, its strnctnre near 
^[ 2 (Af-m)] _ jg same as its strnctnre near = A, which corresponds to a = 0. 

Expanding aronnd —A gives Eq. ([C^) . 
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